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p ■ Abstract 

JL , We match the known chiral perturbation theory representation of the 

qj I TTTT scattering amplitude to two loops with a phenomenological description 

i-^ ' that relies on the Roy equations. On this basis, the corrections to Wein- 

berg's low energy theorems for the S'-wave scattering lengths are worked 

* ^^ I 

k> , out to second order in the expansion in powers of the quark masses. The 

'1-1 ; resuhing predictions, ag = 0.220 ± 0.005, ag = -0.0444 ± 0.0010, contain 

remarkably small uncertainties and thus allow a very sensitive experimen- 
tal test of the hypothesis that the quark condensate is the leading order 
parameter of the spontaneously broken chiral symmetry. 
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1. Introduction 

In the chiral limit, where the masses of the two lightest quarks are turned off, 
QCD acquires an exact SU(2) xSU(2) symmetry. We rely on the standard picture, 
where it is assumed that this symmetry is broken spontaneously and that the 
quark condensate (0|fZM|0) represents the leading order parameter. The quark 
masses act as symmetry breaking parameters, which equip the Goldstone bosons 
- the pions - with a mas^: M^ = 2mi? + 0{m?'), where B stands for the value 
of KOImm |0)|/F^ in the chiral limit. 

Since Goldstone bosons can interact only if they carry momentum, the vrvr 
S'-wave scattering lengths vanish in the chiral limit. Hence these quantities rep- 
resent a sensitive probe of the symmetry breaking generated by the quark masses. 
Weinberg's low energy theorems fl]] state that their values are related to the pion 
mass, which also represents a symmetry breaking effect: 

7 A/f 2 A/f2 

° 327rF2 ^ V ; , i67rF2 ^ ^ > ^ > 

The relations hold to leading order in the expansion in powers of m. The correc- 
tions of order m? were worked out in 0. As an example, we quote the result for 
the combination 2a^ — Soq, which is particularly simple: 

2«o - 5a^ = ^ {l + \Ml{r^), + %^} + O(m^) . (2) 

n 

{r'^)s is the mean square radius of the scalar form factor and ^ measures the pion 
mass in units of the scale set by AirF.^ ~ 1 GeV: 

With the evaluation of the chiral perturbation series to two loops described in 
ref. 10, the low energy expansion of the scattering amplitude is now known to 
next-to-next-to-leading order. The purpose of the present paper is to analyze the 
consequences for the scattering lengths. 

2. Chiral representation of the scattering amplitude 

The two-loop representation yields the first three terms in the low energy expan- 
sion of the partial waves: 

tlis) = t\{s), + t\{s), + t\{s), + 0(/) . (4) 



^We disregard isospin breaking and set rriu = rud = m. In the numerical work, we identify 
Mtt with the mass of the charged pion and use F^r — 92.4 MeV. 



Since inelastic reactions start showing up only at O(p^), unitarity implies 

Im4(.)=a(.)|4(.)P + 0(/) , a(.) = (l-4M,V^)^ . (5) 

The condition fixes the imaginary parts of the two-loop amplitude in terms of 
the one-loop representation. At leading order, the scattering amplitude is linear 
in the Mandelstam variables, so that only the S- and P-waves are different from 
zero. Unitarity therefore implies that, up to and including O(p^), only these 
partial waves develop an imaginary part. Accordingly, the chiral representation 
of the scattering amplitude can be written in the form 0] 

A{s, t, u) = C{s, t, u) + 3277 {I [/°(s) + I (s - m) U\t) + | (s - t) U\u) 

+i U\t) + \ U\u) - \ U\s)] + 0(p«) , (6) 

where C(s, t, u) is a crossing symmetric polynomial, 

C(s, t, m) = Ci+ SC2 + S^ C3 + (t - uf C4 + S^ C5 + S (t - uf Cg . (7) 

The functions U^{s), U^{s) and f/^(s) describe the "unitarity corrections" asso- 
ciated with s-channel isospin / = 0, 1,2, respectively. In view of lin.t\{s)Q oc s^, 
several subtractions are needed for the dispersive representation of these functions 
to converge. We subtract at s = and set 

, _ ^ r , , a{s')t\s'),{t\s'), + 2Ret\s'),} 

The subtraction constants are collected in C{s,t,u). As only the S- and P- 
waves enter, we have dropped the lower index, {t°,t^,t^} = {tQ,t\,tl}. For 
kinematic reasons, the integrand of the P-wave differs from the one of the 5- 
waves: {eo, ei, €2} = {0, 1, 0}. It is straightforward to check that the result of the 
two-loop calculation p is indeed of this structure. 

3. Effective coupling constants 

In the following, the two- loop result for the polynomial part of the amplitude plays 
a key role. It involves the coupling constants occurring in the derivative expansion 
of the effective Lagrangian, C^ff = C2 + C4 + CQ + . . . The corresponding formulae, 
which specify how the coefficients Ci, . . . , Cq depend on the quark masses, can 
be worked out from the two-loop representation given in ref. (for explicit 
expressions, see 0]). These formulae, in particular contain Weinberg's low energy 
theorems, which in this language state that the leading terms in the expansion 
of the first two coefficients are fixed by M^ and F.,,: ci = —M^/F^ + . . . , C2 = 
l/F^ + ... At first order, the constants ^1,^2,^35^4 from £4 enter, and at 
second order, the chiral representation of the scattering amplitude involves the 
couplings ri, . . . , rg from Cq [^ . We need to distinguish two categories of coupling 
constants: 



a. Terms that survive in the chiral limit. Four of the coupling constants that 

enter the two-loop representation of the scattering amplitude belong to this 
category: ii, £2, r^, tq. 

b. Symmetry breaking terms. The corresponding vertices are proportional to a 

power of the quark mass and involve the coupling constants £3, £4, ri, r2, 

The constants of the first category show up in the momentum dependence of the 
scattering amplitude, so that these couplings may be determined phenomenolog- 
ically. The symmetry breaking terms, on the other hand, specify the dependence 
of the amplitude on the quark masses. Since these cannot be varied experimen- 
tally, information concerning the second category of coupling constants can only 
be obtained from sources other than vrvr scattering. The constants r„ from £5 
only generate tiny effects, so that crude theoretical estimates suffice, but the cou- 
plings £3 and £4 from £4 do play an important role and we now briefly discuss 
the information that we will be using for these. 

The crucial parameter that distinguishes the standard framework from the one 
proposed in ref. ||^ is £3. This coupling constant determines the first order correc- 
tion in the Gell-Mann-Oakes-Renner-relation: M^ = 25m {l - ^^£3 + 0(^2) |_ 
The value of £3 is not known accurately. Numerically, however, a significant 
change in the prediction for the scattering lengths can only arise if the crude 
estimate in ref. 0, 

£3 = 2.9 ±2.4, (9) 

should turn out to be entirely wrong. We do not make an attempt at reducing 
the uncertainty in £3 within the standard framework, but will explicitly indicate 
the sensitivity to this coupling constant. 

Chiral symmetry implies that the coupling constant £4 also shows up in the 
expansion of the scalar radius in powers of the quark masses [0] : 

(-')^ = ^{^"4-§+eA.+o(a} • (10) 

As pointed out in ref. , the scalar radius can be determined on the basis of a 
dispersive evaluation of the scalar form factor. We have repeated that calculation 
with the information about the phase shift Sq{s) obtained in ref. [0. In view of 
the strong final state interaction in the S-wave, the scalar radius is significantly 



larger than the electromagnetic one, (r^)^^ = 0.439 ± 0.008 fm^ [|TT|. The result, 

(r^), = 0.61 ±0.04 fm^ , (11) 

confirms the value given in ref. |0 and is consistent with earlier estimates of 
the low energy constant £4 based on the symmetry breaking seen in Fk/F.^ or 
on the decay K -^ Triu |T2[, but is considerably more accurate. Since the chiral 



representation of the scalar form factor is known to two loops [0 , the dependence 
of the correction A^ on the quark masses is also known. In addition to ii, . . . , £4, 
the explicit expression involves a further term, rs2, from £5 El- I^ ^^^ following, 
we use this representation to eliminate the parameter £4 in favour of the scalar 
radius. 

4. Low energy theorems 

We now show that chiral symmetry implies two relations among the coefficients 
Ci, . . . , C4. For this purpose we consider the combinations 

C,^F^{c2 + 4M'Ac3-c,)} , C2^^{-c^+AMt{c,-c,)}. (12) 

n 

Chiral symmetry implies that, if the quark masses are turned off, both Ci and C2 
tend to 1. The contributions from C3 and C4 ensure that the first order corrections 
only involve the symmetry breaking couplings £3 and £4. Eliminating £4 in favour 
of the scalar radius, the low energy theorems take the form 

C^i = l + ^(r^L + ||+e^A, + 0(a , (13) 



At first nonleading order, Ci is fully determined by the contribution from the 
scalar radius, while C2 also contains a contribution from £3. Inserting the values 
{r'^)s = 0.61 fm and £3 = 2.9 and ignoring the two-loop corrections Ai, A2, 
we obtain Ci = 1.103, C2 = 1.117. The value of C2 differs little from Ci - as 
stated above, the estimate (|^) implies that the contributions from £3 are very 
small. The size of the two-loop corrections will be discussed later on, when the 
phenomenological information about the coefficients Ci,... ,04 has been sorted 
out. 

5. Phenomenological representation of the scattering amplitude 

Our predictions for the scattering lengths are based on the comparison of the chi- 
ral representation with the one that follows from analyticity and unitarity alone. 



As shown by Roy ||Tj], the fixed-t dispersion relations can be written in such a 
form that they express the tttt scattering amplitude in terms of the imaginary 
parts in the physical region of the s-channel. The resulting representation for 
^4(5, t, u) contains two subtraction constants, which may be identified with the 
scattering lengths Og and a^. Unitarity converts this representation into a set of 
coupled integral equations, which we recently examined in detail |TD[. The upshot 
of that analysis is that Oq and Oq are the essential low energy parameters: Once 
these are known, the available experimental data determine the behaviour of the 



vrvr scattering amplitude at low energies to within remarkably small uncertain- 
ties. In the present context, the main result of interest is that the representation 
allows us to determine the imaginary parts of the partial waves in terms of Og 
and Oo- Since the resulting representation is based on the available experimental 
information, we refer to it as the phenomenological representation. 

The branch cut generated by the imaginary parts of the partial waves with 
£ > 2 starts manifesting itself only at 0{p^). Accordingly, we may expand the 
corresponding contributions to the dispersion integrals into a Taylor series of the 
momenta. The singularities due to the imaginary parts of the S- and P-waves, 
on the other hand, show up already at 0{p^) - these cannot be replaced by a 
polynomial. The corresponding contributions to the amplitude are of the same 
structure as the unitarity corrections and also involve three functions of a single 
variable. We subtract the relevant dispersion integrals in the same manner as for 
the chiral representation: 

Since all other contributions can be replaced by a polynomial, the phenomeno- 
logical amplitude takes the form 

47r — 

A{s, t, u) = 167fal + ^^ (2a° - Sa^) s + P{s, t, u) (15) 



+327r {lW\s) + |(s - u)W\t) + |(s - t)W\ 



u 



We have explicitly displayed the contributions from the subtraction constants Og 
and Oq. The term P[s,t,u) is a crossing symmetric polynomial 

P{s,t,u)=pi+p2S +P3s'^ +P4{t - U^ +P5S^ +PQs{t - U^ . (16) 

Its coefficients can be expressed in terms of integrals over the imaginary parts 
of the partial waves. We do not list the explicit expressions here, but refer to 
0. In the following, the essential point is that the coefficients pi, . . . ,Pq can be 
determined phenomenologically. 

6. Matching the two representations 

We now show that, in their common domain of validity, the two representations of 
the scattering amplitude specified above agree, provided the parameters occurring 



therein are properly matched. The key observation is that, in the integrals (14), 
only the region where s' is of order p"^ matters for the comparison of the two 
representations. The remainder generates contributions to the amplitude that are 
most of order p^. Moreover, for small values of s', the quantities Imt^(s') are given 

6 



by the one-loop representation, except for contributions that again only manifest 
themselves at 0{p^). This implies that the differences between the functions 
W^{s) and U^ (s) are beyond the accuracy of the two-loop representation ^. 
Hence the two descriptions agree if and only if the polynomial parts do, 

477- _ 

C(s, t, u) = IGnal + ^^ (2a° - Sa^) s + P(s, t, u) + 0(p«) . (17) 

Since the main uncertainties in the coefficients of the polynomial P{s,t,u) arise 
from their sensitivity to the scattering lengths a^, a^, the above relations es- 
sentially determine the coefficients ci, . . . , cg in terms of these two observables. 



The same then also holds for the quantities Ci, C2 defined in eq. ([T^). The 
corresponding low energy theorems for a^ and Oq are of the form 

4 = -^^ + Mt ao + 0{m') , 4 = -^ + M^ a, + O^) , (18) 

with Cq = y(12 Ci — 5 C2). The terms ao, «2 stand for integrals over the imaginary 
parts of the partial waves that can be worked out from the available experimental 
information. Since the behaviour of the imaginary parts near threshold is sensitive 
to the scattering lengths we are looking for, the same applies to these integrals. 
In the narrow range of interest, the dependence is very well described by 

M>o = -04478 + .30Aa° - .37Aal + .5(Aa°)2 - 1.2Aa°Aa2 + 1.8(Aa^)2 , 
M>2 = .0055 + .023Aa° - .095Aa^ + .01(AaJ^)2 - .12AaJ^Aa^ + M{Aalf , 

with Aa[! = a^o- 0.225, Aag = ^0 + 0.03706. 

7. Results for Oq and Oq 

The representation (|1^) splits the correction to Weinberg's leading order formu- 
lae into two parts: a correction factor C„, which at first nonleading order only 
involves the scalar radius and the coupling constant £3, and a term a„ that can 
be determined on phenomenological grounds. 

Inserting the one-loop prediction for Ci, C2 in the relations (|18|) and solving 
for Og, Oq, we obtain the following first order results: 

a° = 0.2195, 0^ = -0.0446, 20^-50^ = 0.662. (19) 

The two-loop corrections Ai and A2 involve the coupling constants £1, £2, ^3, 
the scalar radius, as well as the terms ri, . . . , r4, rs2 from £5. The size of the con- 
tributions from the latter may be estimated with the resonance model described 
in refs. [0, |T3|. The constants £1, £2 can then be determined numerically with the 
phenomenological values of C3 and C4. The resulting two- loop corrections for the 
scattering lengths are very small. The numerical result is sensitive to the value 



of the scale /i at which the renormahzed couphng constants r^(/i) are assumed 
to be saturated by the resonance contributions. For 500 MeV < ^ < 1 GeV, the 
corrections vary in the range -0.001 < ^^Ai < 0.003, -0.004 < ^^As < 0.001. 
In the following, we use the resonance model at the scale /x = Mp and take the 
above range as an estimate for the uncertainties to be attached to the two-loop 
corrections. 

To estimate the errors due to the phenomenological input, we vary the imag- 
inary parts within the range discussed in ref. [1^, where the Roy equations are 
used to determine the behaviour of the S- and P-waves below 800 MeV. If the 
S'-wave scattering lengths are held fixed, the variations in the values of ao, ^2 
are dominated by those in the input used for the phases at 800 MeV. The corre- 
sponding uncertainty in {a^, Oq, 2aQ — Soq} amounts to {±.9, ±.2, ±1} • lO"'^. 

We conclude that the uncertainties are dominated by those in £3 and (r^)^. 
Adding the remaining sources of error up, we obtain^ 

a[] = 0.220 ± 0.001 - 0.0017 A^ + 0.027 A^2 , (20) 

al = -0.0444 ± 0.0003 - 0.0004 Afg - 0.004 A^2 , 

with I3 = 2.9 + A4, (r^)^ = 0.61fm2(l + A^2). Inserting the estimates (|), (|n|), 
we arrive at our final result: 

ag = 0.220 ± 0.005 , ag = -0.0444 ± 0.0010 , . . 

2a[j - 5ag = 0.663 ± 0.006 , ag - ag = 0.265 ± 0.004 . ^ ^ 

8. Discussion 

The truncation of the chiral perturbation series represents an inherent limitation 
of our calculation. We have shown, however, that the corrections of 0{p^) barely 
change the predictions for Og and Cq obtained at 0{p'^). For this reason, we expect 
the contributions from yet higher orders to be entirely negligible. 

The rapid convergence of the series is a virtue of the specific method used 
to match the chiral and phenomenological representations. To demonstrate this, 
we briefly discuss the alternative approach used in ref. [^ |^, where the results 
for the various scattering lengths and effective ranges are obtained by directly 
evaluating the chiral representation of the scattering amplitude at threshold. We 
instead express the amplitude in terms of three functions of a single variable s 
and match the coefficients of the Taylor expansion at s = - in this language, the 
approach of ref. 0, |^ amounts to a matching at threshold. It is straightforward 
to work out the chiral representation at threshold with the values of the effective 
coupling constants that we find with our method. Truncating the series at O(p^), 
we obtain {a[],ao,2a[] - 5al} = {0.200,-0.0445,0.624}, in agreement with the 
2 



The tiny errors given in eg. (|20D merely indicate the noise seen in our calculation - we do 



not claim to describe the scattering amplitude to that accuracy (compare section 14.1 of 10 ). 



valuesQ of ref. 0: {0.20 ± 0.01, -0.042 ± 0.002, 0.614 ± 0.018}. At threshold, 
the terms of 0{p^) are by no means neghgible: They take the values obtained at 
0{p'^) into {0.215, —0.0445, 0.653}, in agreement with the results of ref. [^ |15], 0. 
These numbers describe the expansion of the scattering lengths in powers of the 
quark massesQ to order ttv'. In the case of Oq, for instance, the numerical values 
to order m, m? and m? are Oq = 0.159, 0.200 and 0.215, respectively - these 
correspond to the diamonds shown in the figure. 

The reason why the straightforward expansion of the scattering lengths in 
powers of the quark masses converges rather slowly is that these represent the 
values of the amplitude at threshold, that is at the place where the branch cut 
required by unitarity starts. The truncated chiral representation does not de- 
scribe that singularity well enough, particularly at one loop, where the relevant 
imaginary parts stem from the tree level approximation. The matching must be 
done in such a manner that the higher order effects are small. In contrast to 
a matching at threshold - that is, to the straightforward expansion of the scat- 
tering lengths - our method fulfills this criterion remarkably well: We are using 
the expansion in powers of the quark masses only for the coefficients Ci and C2, 
while the curvature generated by the unitarity cut is evaluated phenomenologi- 
cally. Solving eq.(|18D with the expansion of these coefficients to order 1, m and 
m^, we obtain a much more rapidly convergent series: Og = 0.197, 0.2195, 0.220. 

If the effective coupling constants are the same, the only difference between the 
two approaches is the one between the functions W^{s) and U\s). In particular, 
the results for Oq, Oq only differ because the numerical values of W^s) and U\s) 
at s = 4M^ are not the same. As mentioned above, the difference between the 
two sets of functions affects the scattering amplitude only at O(p^) and beyond. 
Numerically, however, it is not irrelevant which one of the two is used to describe 
the effects generated by the unitarity cuts: While the functions W\s) account 
for the imaginary parts of the 5*- and P-waves to the accuracy to which these 
are known, the quantities U\s) represent a comparatively crude approximation, 
obtained by evaluating the imaginary parts with the one-loop representation. 

9. Conclusion 

Our analysis relies on two ingredients: the evaluation of the chiral perturbation 
series for the vrvr scattering amplitude to two loops ^ and the phenomenological 
representation obtained in ref. |jlO[ by solving the Roy equations. We have shown 



■^As stated in g], the error bars given only measure the accuracy to which the first order 
corrections can be calculated; they do not include an estimate of contributions due to higher 
order terms. The small numerical differences arise partly from the manner in which the coupling 
constants £1, £2 are determined, partly from the values used for F^^ - the old number, 93.3 MeV, 
does not account for radiative corrections. 

"^More precisely, the expansion parameter is the physical pion mass and, in lieu of the 
coupling constant F, the physical decay constant is held fixed. 



that the comparison of the two descriptions allows us to predict the S'-wave 
scattering lengths at the 2-3% level of accuracy. 




Scattering lengths: theory versus experiment. The shaded region represents the inter- 



section of the domains allowed by the old data and by the Olsson sum rule, see |]iO| 
for details. The ellipse indicates the impact of the new, preliminary K^^ data [18|. 
The cross shows the result of the present paper. The three diamonds illustrate the 
convergence of the chiral perturbation series at threshold (see text). The one at the 
left corresponds to Weinberg's leading order formulae. The two black solid lines are 
the boundaries of the universal band llC 



We emphasize that our result (21) relies on the standard picture, according to which 
the quark condensate represents the leading order parameter of the spontaneously 
broken symmetry. The scenario investigated in ref. ^ concerns the possibility that 
the Gell-Mann-Oakes Renner formula fails, the second term in the expansion M^ = 
2Bm<l — 2C^3 + 0{^'^) \ being of the same numerical order of magnitude or even 
larger than the first. Note that for this to happen, the value of \l^\ must exceed the 
estimate (|9|) by more than an order of magnitude. 



The constraints imposed on Og, 



On by the available experimental information are 



shown in the figure. The ellipse represents the 68 
by combining the new, preliminary K^.^ data ||l 



^ confidence level contour obtained 
with earlier experimental results. 
Concerning the value of Cq, the ellipse corresponds to the range 0.2 < Oq < 0.25. The 
representation (^) shows that this range only puts very weak limits on the value of 
I^. The precision data on the reaction -kN — >■ TrvriV near threshold [l^ provide an 
independent measurement of Og. Unfortunately, however, the systematic errors of the 



10 



Chew-Low extrapolation that underhes this determination do not appear to be under 



good control [gO[ - for a detailed discussion, we refer to ||21| 



The figure shows that the values of Og and a^ are strongly correlated. The cor- 
relation also manifests itself in the Olsson sum rule |l^, which according to ref. |10| 
leads to 2a^ — Sag = 0.663 ± 0.021 + l.lSAag — 1.01Aa57in perfect agreement with our 
result in eq. (]2l|). Note, however, that this combination is not sensitive to £3 - accurate 
experimental information in the threshold region is needed to perform a thorough test 
of the theoretical framework that underlies our calculation. We are confident that the 



forthcoming results from Brookhaven [|^, CERN |22, 23] and Frascati [^] will provide 
such a test. 
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